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Abstract 

The extraordinary ability of space-charge waves in plasmas to accelerate charged particles at gra- 
dients that are orders of magnitude greater than in current accelerators has been well documented. 
We develop a phenomenological framework for Laser WakeField Acceleration (LWFA) in the 3D 
nonlinear regime, in which the plasma electrons are expelled by the radiation pressure of a short 
pulse laser, leading to nearly complete blowout. Our theory provides a recipe for designing a LWFA 
for given laser and plasma parameters and estimates the number and the energy of the accelerated 
electrons whether self-injected or externally injected. These formulas apply for self-guided as well 
as externally guided pulses (e.g. by plasma channels). We demonstrate our results by presenting 
a sample Particle-In-Cell (PIC) simulation of a 30/sec, 200TW laser interacting with a 0.75cm 
long plasma with density 1.5 x 10 18 cm -3 to produce an ultra-short (10 fs) mono-energetic bunch 
of self-injected electrons at 1.5 GeV with 0.3nC of charge. For future higher-energy accelerator 
applications we propose a parameter space, that is distinct from that described by Gordienko and 
Pukhov [Physics of Plasmas 12, 043109 (2005)] in that it involves lower densities and wider spot 
sizes while keeping the intensity relatively constant. We find that this helps increase the output 
electron beam energy while keeping the efficiency high. 

PACS numbers: 
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In plasma based acceleration a laser or particle beam creates a plasma wave wakefield 
with a phase velocity close to the speed of light, c. The acceleration gradients in 

these wakefields can easily exceed 10 GeV/m which is nearly three orders of magnitude 
larger than that achieved in conventional RF technology. A particle injected in such a wave 
with sufficient initial energy can interact with the longitudinal component of the electric 
field for a long enough time, that its energy gain is significant. 

An intense laser pulse or electron beam can expel the plasma electrons outward to create 

i n 

a bare ion column |3j, |4| which has ideal accelerating and focusing properties for electron 
beams Q] . A recent important development in this field has been the simulation and exper- 
imental observations of these highly-nonlinear 3D wakefield structures and the acceleration 
of e.ectrons fl fl B, fl fl Q . In the laser-driven cases the y we r e shown to ~ q nasi- 
monoenergetic beams of electrons with energies on the order of 100 MeV QflQ- In these 
experiments, the energetic electron generation mechanism involves an interplay of complex 
phenomena, namely: the evolution of the laser, the self-trapping of electrons and their 
influence on the wake. 

In this article, we describe how to extend the results from the recent experiments towards 
a more stable regime conducive to making an accelerator. We will show that by using near 
term lasers in the 0.1 to 3PW range it will be possible to generate of 1 — 13GeUmono- 
energetic electron beams with nC of charge in a single stage without the need for external 
guiding. Using external guiding energies up to 120GeV can be achieved in a single stage for 
a 3PW laser. Unlike in the recent experiments, this regime is characterized by lasers whose 
intensity and spot size are matched and by relatively low plasma densities. The matched laser 
profile (both transverse and longitudiunal) evolves little during the acceleration distance 
which greatly exceeds the Rayleigh length. 

We are motivated to determine a path towards making a compact electron accelerator 
based on the LWFA concept for use in high energy physics or as a light source. Efficiency and 
beam quality are key issues one should keep in mind when making such designs. As shown 
by Katsouleas et al. H[ , in order to achieve sufficient beam loading efficiency (from wake to 
particles) and to reduce the beam emittance growth, the spot size of the wake (hence laser) 
should be less than a plasma skindepth. In addition, in order to most efficiently transfer 
the laser energy to the wake over the entire effective accelerating distance (pump depletion 
length matching dephasing length), the laser normalized vector potential (laser intensity) 
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must be larger than unity. These two conditions are based on linear theory, however, taken 
together they imply that LWFA needs to be operated in a 3D nonlinear regime, in which 
the laser power exceeds the critical power for relativistic self-focusingj^]. 

It was shown by Sun et al. [3] that for intense and narrow lasers the transverse pondero- 
motive force can lead to complete electron blowout (cavitation) when the laser power slightly 
exceeds the critical power for relativistic self-focusing,P c , and the laser is focused tightly. 
This analysis was done for long pulses (although it assumed the ions were fixed). Mora 



and Antonsen |l2j showed in reduced PIC simulation that cavitation could occur for power 
above P c . When using ultra-short pulses, cr < wo — 2^ c/u p , where ao = eA^ aser j \mc 2 ) 
for a linearly polarized laser and wq the laser spot size, the ions are indeed stationary, and 
a region nearly void of electrons (ion channel) is generated behind the laser pulse. This 
channel exerts an attractive Coulomb force on the blown-out electrons causing them to rush 
back toward the axis thereby exciting a wakefield. In 1991 Rosenzweig et al. |4J showed that 
wakes excited in this nonlinear blowout regime (they used a particle beam driver) had ideal 
focusing and accelerating properties. Very soon thereafter it was shown using 2D simulations 
that similar wakes could be made by lasers 

In 2002 Puhkov and Meyer-ter-vehn used 3D particle-in-cell simulations of LWFA at very 
high laser intensities and plasma densities without external guiding and found that electrons 
were self-injected in the blowout regime and that they formed a quasi mono-energetic beam 
as they were accelerated [5(] . In their simulations the blowout cavity formed a sphere leading 
them to subsequently call this the "bubble" regime for LWFA 141 ] . Recently, a theoretical 
analysis 15|, 116J showed that for sufficiently intense and ultra-short lasers, ao > 4, the 
resulting electron density structure will resemble a spherical cavity and it was shown that 
the laser spot size needed to be matched to the maximum blowout (bubble) radius in order 
to achieve the most ideal wakes. It is worth noting for 2 < a < 4 electron blowout still 
occurs with the cavity slightly deviating from a spherical shape. 

The recent experimental results are promising but they still rely on substantial evolution 
of the laser and the energy spread and emittances are not sufficient for high-energy physics 
applications. What paths are available for extending these experiments towards higher 
energies with more control using the "bubble" regime? The original simulation results on the 
"bubble" regime in ref. |5( were based on a simulation of a 3607W laser propagating through 
a 0.01 critical density (1.74 x 10 19 cm -3 for A = 0.8«m) plasma. Do the key features of this 
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simulation remain for different laser powers, pulse lengths, spot sizes, and plasma densities? 
Answering these questions is difficult because the inherent physics is highly nonlinear. 

One approach was put forth by Gordienko and Puhkov (GP) 17]]. They argued that 
for a >> 1 the speed of all plasma electrons is very close to the speed of light and under 
this condition all quantities will scale with a single similarity parameter, S = to 
some power. The coefficients in front of the scalings are determined from simulations (or 
experiments). The resulting expressions are therefore only strictly valid so long as the 
laser's transverse and longitudinal profile, aspect ratio etc. remain the same. Another 
approach presented here is to use a phenomenological description. We identify the important 
physics as wake excitation (amplitude and phase velocity), pump depletion (pulse evolution), 
dephasing (between particles and wake), and beam loading. We use these concepts to develop 
expressions for predicting the number of electrons, the electron energy, and overall efficiency. 

For any useful accelerator, the accelerating structure's (wake) amplitude, wavelength, and 
phase velocity need to be reasonably stable during the entire acceleration process. We show 
that a stable wake can be excited in the blowout regime when the laser's spot size is roughly 
matched to the blowout radius. We also show that under these conditions expressions for 
the wake amplitude and phase velocity (dephasing length), the laser pump depletion length, 
and the number of accelerated electrons can then be obtained including their coefficients 
based on simple physics arguments and previously published ideas. Importantly, unlike 
when using a similarity theory, the coefficients derived from the phenomenological theory 
can be extended to propagation in plasma channels as well to when an external beam of 
electrons is injected. 

We will show that while there is overlap in the ideas put forth here and by GP there are 
also profound differences besides those just described. Their analysis can only be valid when 



a o ^ 2-y/n c /n p while ours is for a > 2. Their supporting simulations used a = 20 — 80 
(for each direction of a circularly polarized laser which is the same as 28 — 113 for a linearly 
polarized laser) and n/n c = 0.02 — 0.08. Their analysis appears to ignore dephasing leading 
to different scaling laws. We will show that in contrast to the work of GP, our expressions 
advocate using lower ao and lower densities. In fact, we will argue that in order to obtain 
controlled acceleration one should operate with an a Q too low for the key assumption in the 
theory of GP to even apply. Despite these differences we borrow the term "bubble" from 
ref. 14j to describe the near spherical ion channel behind the laser pulse. 
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To verify our results, we have carried out many computer experiments using the 3D, 



particle-in-cell code OSIRIS 



18] , that explore a wide range of plasma densities, laser powers, 



and spot sizes (see later). We highlight one simulation that is very relevant to near term 
experiments. Details of the simulation parameters and setup are provided in the Appendix. 
In this simulation, a 30 fs (FWHM) 0.8/zm laser pulse containing 200 TW of power is focused 
to a spot size wq = 19.5/xm at the entrance of a 1.5 x 10 18 cm -3 density plasma to give a 
normalized vector potential of ao = 4. The laser is circularly polarized (with normalized 
vector potential 4/ a/2 in each direction) and has a Gaussian transverse profile. The plasma 
is 0.75cm long which corresponds to more than 5 Rayleigh lengths. 

We next develop a phenomeno logical theory for LWFA in the blowout regime. A light 
pulse can be guided with a nearly constant (matched) spot size by a plasma channel that has 
a parabolic refractive index/ density profile with a maximum/minimum on axis. The index of 
refraction in a plasma can be expanded asQ, 3 V = ck/u ~ 1 — ^1 + ^^fjz + ^ ~ -§■ 
where An c parametrizes an external density channel, An is a density depletion from the 
transverse ponderomotive force and the term is due to relativistic mass corrections. 
The characteristic density change required to optically guide such a profile with little spot 
size oscillation is An c = l/(7rr e Wg) 19j, where r e = e 2 /mc 2 is the classical electron radius 
and wo is the laser spot size. If the density depression is normalized to the plasma density 
n p , this condition becomes, An c /n p ~ 4/(k p Wo) 2 . An equivalent change to the index of 
refraction from relativistic mass corrections can also self-guide a laser if a^/8 > 4/(/c p tf ) 2 
or P > P c , where P c = YI uj 2 / uj 2 [GW] is the critical power for relativistic self focusing 

Unfortunately, even when P > P c and the spotsize is sufficiently small, the index of 
refraction needs a distance on the order of c/u p to build up. Therefore, it is often thought 
that a short pulse laser r < l/u) p cannot be self-guided and some form of external optical 
guiding is needed [3]. However, as described in ref. 21] for P/P c » 1 a degree of self- 
guiding for short pulses is possible because the leading edge of the laser locally pump depletes 
before it diffracts and the back of the pulse is still guided in the ion column region. In case 
of external optical guiding P/P c » 1 is not necessary, instead P/P c > 1 can be used which 
also leads to clear blowout. 

The stable self-guiding of an intense short pulse is illustrated in figure [H We find that for 
self-guided propagation of the laser beam, without significant variations of the pulse profile 
over the interaction distance, its spot size and intensity must be appropriately chosen. A 
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rigorous derivation for the matched spot size and the laser profile is not currently available. 
Still, as noted in ref. , [ijj] the requirements for a matched profile can be estimated by 
assuming that the transverse ponderomotive force of the laser /cpVajj/7 ~ ao/(k p R), where 
R is the blowout radius, is roughly balanced by the force of the ion channel E r ~ k p R, which 
pulls back the ponderomotively expelled electrons. Equating these two expressions yields 
k p R ~ y/cbQ. Furthermore, the blowout radius scales with the laser spot size, k p wo ~ k p R. 
Through simulations we have found that a more refined condition which leads to only slight 
oscillations in the spot size is: 

k p R ~ kpWo = 2t/oo (1) 
We find that this relationship still holds for a > 2. In the work of Gordienko and Pukhov 



171 ] it was assumed k p R ~ k p wo = 1.12^/ao (we rewrote this for a linearly polarized laser). 



We can use Eq. [T] to reformulate the matched beam spot size condition as: 

a ~ 2(P/P C ) 1 / 3 (2) 



As alluded to earlier, the laser etches back due to local pump depletion. In ref. [21(, an 
estimate of the etching rate, i.e., the etching velocity, based on nonlinear ID effects was 
given as v etc h — cuj p /u}q] and it was verified in 2D PIC simulations. We find that on average 
this estimate also agrees very well with the observations in numerous 3D PIC simulations. 
Therefore, the laser will be depleted after a distance (pump depletion length) 

c 

Letch — CT FWHM — -^CT FWHM (3) 

Vetch 

The front of the laser, that excites the wake, moves backward as the pulse etches back with 
Vetch- The phase velocity of the wake can therefore be expressed as ~ v g — v etc h, where 
v g is the linear group velocity of light in a very underdense plasma Uq; therefore 

v,/, ~ c[l — 3UL 1 / (2o>q)]. The distance that the trapped electrons travel until they outrun the 
wave (dephasing length) is: 

R^l^LR (4) 



c — 3 cj. 



We find in numerous ID, 2D and 3D simulations that the etching velocity and hence this 
dephasing estimate works well for 2 < ao < 2^J^. The estimate for the upper value of a® 



is discussed later. 
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For the energy gain we may write the obvious equation 

AE = qE LW L acc = e LW l acc mc 2 (5) 

where Elw is the average accelerating field of the beam loaded wake, L acc is the acceleration 
length, €lw = gElw / \mao p ) and l acc = ui p L acc /c. The desired acceleration length is the 
dephasing length, so we impose the condition L etch > L d => ct fwhm > 2R/3. If the pulse is 
too short dephasing will not be reached, and the electron beam may have significant energy 
spread. We note that this condition is approximate, since the laser starts diffracting as 
soon as its intensity is insufficient to sustain self focusing and this happens before the pulse 
is completely pump depleted. Additionally the injected particles need to slightly pass the 
dephasing point (phase space rotation) so that the energy spread is minimum. On the other 
hand, the length of the pulse should not be too large, because the laser field could interact 
with the trapped electrons and degrade the beam quality. 

We need an expression for eL,w which is the average accelerating field experienced by an 
electron, (E z ). In ref. [ljj] it was shown that for a > 4 the ion column formed a sphere 
and that the accelerating field, for the most part, depends linearly on the distance from the 
middle of the sphere (figure EJa)). This is confirmed in figures l(a)-(c) where a bubble void 
of electrons roughly forms a circle and in figures 2(a), (b) where a lineout of eE z /(mcuj p ) 
along the axis is shown. We also found by theory and 3D PIC simulations that a spherical 
shape is still roughly formed for 2 < a < 4 with matched laser profiles. Interestingly, 
although the physics is very different, the nonlinear wakes formed in ID also exhibit a linear 
slope of the same amount. 

Because the bubble is roughly a sphere and the electrons are either self or externally 
injected at the rear, the electrons then travel a relative distance R before they dephase. The 
peak useful accelerating field is eE z , max / {mcu) p ) = tJEq and because the wakefield is roughly 
linear, the average field is half of the peak: €lw = eE ZtTnax / \2mcuj p ) ~ ^fa^jl . We can 
therefore write the approximate equation for the energy gain: 



S 



AE ~ ^mc 2 (— ^ a 



,2 



1/3 / \ 2/3 



~ mC " 1 2 5/ 2 — 



100 J Wcm- 3 ]J VV/H/ 
We emphasize the much stronger dependence of the beam energy on the plasma density 
than on the input laser power. However, when the plasma density is lowered for fixed power, 
ensuring self-guided propagation of the leading edge of the laser is more challenging. This 
can be accomplished by plasma channels or to some degree by self-guiding. As we argue 
later, for self-guiding to occur P/P c needs to increase as the plasma density decreases. We 
can rewrite equation (EJ) in terms of the critical power for relativistic self-focusing, P c : 

On the other hand, if a plasma channel is used, P/P c can be kept as low as 1. As shown 
before, the channel depth An c /n p needed is 4/(k p w ) 2 ~ l/a . So as long as a > 2 is used, 
the normalized channel depth is small (less than 0.5). It is also worth noting here that when 
a channel is used, the channel parameters (width and depth) should be chosen based on the 
matching condition for given laser power and plasma density (Eqs. [I] and [2]) so that leading 
front of the laser is guided by the density channel and the back of the laser is guided by the 
matched ion channel. 

The electron beam surfing on the wake can be either self-injected as shown in our sample 
simulation (figure QJb)) or externally injected from some other source. For self-injection, 
particles in the rear of the blowout region must be able to catch up with the wake. The 
physical condition for this to happen is twofold: first, the blowout radius should be large 
enough so that when the particles reach the rear of the bubble, they move predominately 
in the forward direction with speed close to the speed of light. Second, at the rear portion 
of the ion channel, trajectory crossing occurs leading to a narrow sheath with the highest 
accelerating and focusing fields. Therefore, even though electrons initially have a 7 (energy) 
substantially below the wake's Lorentz factor 7^, they can easily achieve sufficient energy as 
they are accelerated while they slowly drift backwards (relatively to the pulse) in the sheath. 
In our sample simulation, the effective 7 of the wake is around 20 and the normalized blowout 



radius is around 4. The initial energy 7 of those trapped electrons is substantially smaller 
than 20. For even lower plasma densites, we have performed a number of simulations, where 
an electron beam with 7 exceeding 10000 was used as the driver instead of a laser and we 
observed self-injected electrons in each case for a normalized blowout radius around 5. This 
indicates that for laser wavelengths in the 0.8/im range and plasma densities of interest, 
self-injection will always happen when we keep the normalized blowout radius around 4 ~ 5. 



14[ that v /oq > 70, or ao > 400 for our 



This differs significantly from the claim in ref. 
sample simulation, for self-trapping to occur. 

In the regime presented here the self-injected electron bunches are highly localized in 
space with a half- width of the first bunch of only ~ 10/s, i.e. lc/u p . Once a sufficient 
number of electrons have been trapped the trapping process terminates, as seen in figure 
[T^c). The first electron bunch reaches an energy of 1.5 GeV and its energy spectrum is 
presented in figure [2]^c). The normalized emittances are shown in figure [31 They may be 
estimated as the product of the beam spot size, which roughly scales with 1 / y/n^, with the 
spread in the momentum perpendicular to the acceleration direction, which scales with the 
relativistic ponderomotive potential (~ ao). These simple considerations show that as we 
move to lower densities in order to achieve higher energy particles the emmitances of the 
self-injected electrons will increase. This suggests that for the electron beam to be useful 
for high energy physics or light source, external injection may be more attractive. As an 
interesting aside, simulations also reveal the trapping and acceleration of a second distinct 
bunch in the second bucket (see figure EJ^c)), which has a lower energy because the average 
electric field it experiences is smaller. 

The number, N, of electrons that are accelerated can be estimated from energy balance. 
Hence we examine the partition of field and particle energy within the first bucket. The 
fields inside the ion column have E z , E r , and components. In addition there is kinetic 
energy in the plasma. In a 3D linear or ID nonlinear wake the fields and kinetic energy scale 
together. This is not the case for these 3D nonlinear wakes where an increasing percentage of 
energy ends up in electrons which are blown out well beyond the narrow electron sheath for 
higher laser intensity. Integrating the field energy in the ion channel we find equipartition 
between the energy in the longitudinal field Si and the focusing fields Sf. 
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The trailing particles can recover the field energy in the ion channel and the kinetic 
energy in the narrow electron sheath by changing the shape of the ion channel. The kinetic 
energy in the narrow electron sheath scales the same way as the field energy. By equating 2S 
with the energy absorbed by N particles that travel across the ion channel (we assume the 
average field felt by these particle is E z ^ max /2 and the kinetic energy in the electron sheath 
absorbed by the trailing particles is the same as the field energy), we obtain: 



where a = k p w / (2y/ao) and (3 = k p R/ Using equation ([I]) a ~ 1 ~ (3 we obtain: 



The efficiency scales as the total energy Sf, in the accelerated electron beam (energy gain 
equation ([6]) times particle number from equation (Q) devided by the total laser energy St 
(assuming cr ~ 2y / a c/co , p ): 



which indicates that ao, i.e., (P/P c ) 1//3 cannot be too large if one needs high efficiency. For 
a 200TW, 0.8/zm pulse equation (TTUT) predicts 0.6nC of charge. The charge measured from 
the simulation for the first bunch is 0.3nC. We have also verified these scaling laws by 
monitoring how many electrons can be externally injected before the wake becomes severely 
loaded. 

As shown in figure [T](d), the acceleration process stops before the accelerating bunch de- 
phases. This will not lead to any considerable modifications of the aforementioned formulas, 
particularly because the pump depletion length scales as the dephasing distance and the 
accelerating wakefield decreases as the trapped electrons approach the center of the sphere. 

The beam energy seen in the simulation, 1.5 GeV, is close to that calculated theoretically 
from Eq. ([6]) which is (\6§ =>- AW ~ 1.6GeV. Using formulas ©-(JI]) we see L d ~ 1.31cm > 
0.96cm ~ L p d which is in agreement with our observation that pump depletion happened 
before dephasing. In addition, as the laser pump depleted its rate of diffraction increased 
such that the effective acceleration length was eventually limited by diffraction to a distance 
0.75cm which was less than the pump depletion distance 0.96cm. For a wake with a linear 
slope most of the energy gain occurs near the peak field. This combined with the density 






(10) 



r ~ Sb/Sx ~ i/ a o 



(ii) 



n 



spike at the back of the bucket is why the estimate from eq. ([6]) is still so accurate for this 
non-optimized sample simulation. 

In spite of the complexity of the physics associated with this interaction, the predictions 
by the simple formulas presented in this article are very close to 3D PIC simulation results. 
Good agreement is also achieved between these scaling laws and recent experimental results 



9] despite the fact that the laser powers were slightly below the "threshold" for the 
blowout regime and the laser pulse was not matched transversely nor longitudinally. To be 
in the regime identified by a spherical bubble, one needs do > 4 or equivalently P/P c > 8 and 
cr < which leads to the condition that P > 30(r/30fs)TW. These conditions 

can be relaxed for channel-guided lasers for which P/P c may be smaller. It is then written 
as Oq > 2 or P/P c > 1. We present the comparison between the scaling law for the energy 
Eq. ([6]) and the aforementioned results in figure HI 

The scalings derived above and the underlying physics predict that it is advantageous to 
use moderate intensities and very low plasma densities to increase the output energy and 
keep the efficiency high. The simulations also show that the injection process can be clamped 
and the energy spread of the electron beam is much less in the regime we are proposing, 
which indicates that this regime is also amenable to accelerating externally injected beams 
while maintaining good beam quality. 

In table fl] we compare some of our results to previously published work [IT] , Both 
linear and ID nonlinear formulas fail to describe the 3D nonlinear regime in nearly all 
respects. It is however instructive to compare the scaling laws described above against 



those obtained in ref . [IT 
regime, because in ref. 



using a similarity theory for the high intensity high plasma density 



171 ] the scaling (but not the coefficients) of the fields in the bubble 



and the laser spotsize with ao were identified correctly. The key assumption of ref. 
required that ao >> 1. For the energy gain, they obtained 
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P CT 

Emono ~ 0.65mc 2 W — (12) 

y m A & /e z A 

This formula also implies a dependence on plasma density if cr ~ ^ ~ \^k)c/u p is assumed. 
Under these optimum conditions, formula ( Tl2|) can be rewritten as 

/ D \ 2/3 / \ 1/3 



E mono « 0.16mc 2 — ( : ) ( — ) (13) 
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TABLE I: This table compares formulas and scalings from linear theory, ID nonlinear theory, the 
similarity theory of reference 16 and the 3D phenomenological nonlinear theory in this paper. The 
column labeled ao shows the range of ao that each theory is intended for. k p WQ from linear and 
Id nonlinear theory is kept constant while in both our theory and theory of ref. [13]. k p WQ scales 
with the blowout radius. The acceleration length is limited by dephasing in linear theory and our 
3D nonlinear theory, while it is limited by pump depletion in ID nonlinear theory and the work of 
GP. (Note that in standard ID nonlinear theory the phase velocity of the wake is incorrectly set 
equal to the nonlinear group velocity) Only in our 3D nonlinear theory do dephasing and pump 
depletion scale the same way (if the pulse length is matched to the bubble size). The dephasing 



length is evidently ignored in ref. 13] so the energy gain depends on the pulse length. Linear 
dephasing and Id nonlinear dephasing both scale differently than in the 3d nonlinear regime. The 
pump depletion length from Id nonlinear theory scales as in 3d nonlinear theory but for different 
physics reasons. In ID all the laser energy goes into the wake while in our 3D theory this is not 
the case. We discuss the pump depletion length from ref. l?j] extensively in the text. Both linear 
and Id nonlinear theory incorrectly identify the plasma wavelength and wake phase velocity. The 



17] then the energy scales as 
171 ] all give identical scalings 



resulting energy scalings are different. If one sets uj p t ~ ^/Eq in ref. 
<2q. Therefore linear, ID nonlinear and the similarity theory of ref. 
for the AW with both intensity and density while our results give A scaling as ao to the first 
power. 



We point out the stronger dependence of the energy on the power in these formulas. The 
difference between the scaling laws presented here and those from similarity theory lie in 
the scaling for L acc . We have argued that L acc is limited because the pulse pump depletes 
primarily by giving kinetic energy to the electrons. An electron at the front of the laser is 
pushed forward and to the side. As argued in ref. 2l|] each electron gains an amount of 



13 



energy that scales as (a^/2)mc 2 . Since the laser's energy also scales as this "ID like" 
pump depletion enhances the etching velocity independently of laser intensity. However, for 
extremely large a Q electrons can move forward with a velocity greater than the velocity of the 
leading edge of the pulse; which is at most the linear group velocity because the nonlinearity 
does not develop instantaneously. Equating the forward going velocity of a single electron 
to the linear group velocity gives an estimate for a critical value for the laser amplitude, 



aoc ~ 2^Jn c /n p . We can see that for high plasma density n p /n c « 0.01 ~ 0.08 the critical 
value for a is a 0c ~ 20 ~ 8, but for lower plasma densities n p /n c ~ 0.001, a Qc > 70. 

For laser amplitudes above this critical value, the laser will pump deplete more slowly 
because the kinetic energy given to the each electron will no longer scale as Oq. In this limit a 
constant percentage of the energy may go into the fields of the wake such that T ~ constant 
as predicted by ref. {l7|. The pump depletion length, L p d, can then be estimated from energy 

jrp "Pp op up i up 

balance between the laser and the wakefields, g^-cr ~ -^L pd L pd rs ^f^cr ~ z a oz? CT ~ 
aQL e t c h which does depend on a . Using this scaling for L acc in equation (jSJ) provides the 
same scaling for the particle energy as obtained from similarity theory. However, even for 
this high intensity limit dephasing will still occur because the leading edge cannot move 
faster than the linear group velocity of the highest frequency component, since at the very 
front the amplitude is small 22]. 



Last, it is useful to extrapolate the bubble regime for moderate a to lOGeV and beyond. 
If we keep P/ P c fixed, then the laser power and inverse of the density scale with the desired 
electron energy and the pulse length will scale as the electron energy to the 1/2 power. If 
we scale up our sample simulation, to obtain 15GeV we need a 2PW, 100/s laser and a 
plasma with density of 1.5 x 10 17 cm~ 3 . The only uncertainty is whether self-guiding is still 
possible. 

To obtain condition for self-guiding we resort to quasi-static theory where the index of 



refraction rj ~ 1 - 5 ^ ^ and w , 



l+a? 



- 1 



0. As noted earlier 1 — -p-r ~ ib must 



(i+V) 2 

equal 4/(/c p w ) 2 in order for guiding to occur. For <^ 1, we assume ^ ~ if a2 i therefore 
ip builds up to the necessary value in a distance A£ which scales as A£ ~ (kpWoa ) _1 ~ 
{k p a^ 2 )~ l . The length of the laser that is lost due to diffraction each Rayleigh length (Zr) 
will also scale as (k p a^ 2 )~ l . Self guiding can be achieved if the length lost due to pump 
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depletion in each Rayleigh length also scales the same way: 

Vetch 3 n p A£ j/g 

c z n c Z fi 

This scaling law indicates that as the density is decreased to increase the electron energy 
we need to slowly increase the laser intensity to maintain self-guiding. As a result the 
2PW laser pulse estimate needs to be slightly modified. Instead we choose P/P c — 28 =>- 
ao — 6, and a density n = 4x 10 17 cm -3 . For these parameters stable self guiding should 
occur and lead to acceleration of 1.8nC charge at approximately 9GeV after 10cm of laser 
propagation. Carrying out a full-scale PIC simulation for the 10cm of propagation distance 
is beyond current capabilities. However, we have carried out a 3D PIC simulation for a short 
propagation distance to verify that electrons will still be self-injected. Results are shown in 
figure El 

If both the technology for making meter scale (0.1m ~ lm) low density plasma channels 
for guiding and for synchronized external injection are developed, then much higher energy 
gain can be achieved for the same laser power. For example, using P/P c ~ 1 (and the 
channel depth An c /n p = 0.5), a lOGeV electron beam with ~ 0.6nC charge could be 
obtained with a 250TW, 100/s laser and a 43cm long plasma channel with a minimum 
density n — 1.2 x 10 17 cm -3 . Extrapolating the parameters furrther, we predict that a 
120GeV electron beam with 2.2nC charge generated by using a 3PW, 350/s laser and 18m 
long plasma channel with density n = 1.0 x 10 16 cm~ 3 . 

The work is supported by DOE under grants DE-FC02-01ER41179, De-FG02-03ER54721, 
DE-FG03-92-ER4727 and DE-FG03-NA0065, by NSF under grant Phy-0321345, and by 
FCT (Portugal). The simulations are performed on Dawson cluster, maintained locally by 
UCLA/ATS, and on the IBM SP @ NERSC under mpll3 and gc2. We acknowledge very 
useful conversations with Profs. T. Katsouleas and A. Pukhov. 



I. APPENDIX: PARAMETERS AND SETUP 

The simulation for the 2007W, 30/sec laser pulse in a fully ionized n = 1.5 x 10 18 cm -3 
plasma has a computational window of dimension 101.9 x 127.3 x 127.3/im 3 which moves 
at the speed of light. The number of gridpoints is 4000 x 256 x 256 = 2.62 x 10 8 . The 
resolution in the laser propagation direction z is k Az = 0.2. We assume a preformed fully 
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ionized plasma with uniform density profile. The resolution in the transverse direction is 
kpAx = kpAy = 0.116. We use 2 electrons per cell and a smooth neutralizing immobile 
ion background (the total number of particles is roughly 500 million). A diffraction limited, 
circularly polarized pulse is focused at the plasma entrance with a spot size Wq = 19.5/im 
and the electric field has a symmetric temporal profile of 10r 3 — 15r 4 + 6r 5 where r = 
V2(t — t )/TFWHM- The laser is Gaussian in the transverse direction. The total axial length 
of the plasma is close to 0.75cm, or 300, 000 simulation time-steps. 

The simulation for the 2PW, 100/sec, w = AOfim laser in plasma density n = 4 x 10 17 
has a computational window of dimension 203.8 x 127.3 x 127.3/im 3 which corresponds to 
8000 x 256 x 256 = 5.24 x 10 8 gridpoints with identical resolution as in the aforementioned 
200TW run. Roughly 1 billion particles were used and ran 16000 simulation time-steps. 
The rise-time of the laser was about 2.5 times shorter than the fall-time. 
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FIG. 1: A sequence of 2-dimensional slices (x — z) reveals the evolution of the accelerating structure 
(electron density, blue) and the laser pulse (orange). Each plot is a rectangular of size z = 101.7/xm 
(longitudinal direction, z) and x = 129.3/um (transverse direction, x). A broken white circle is 
superimposed on each plot to show the shape of the blown-out region. When the front of the laser 
has propagated a distance, (a) z = 0.3mm the matched laser pulse has clearly excited a wakefield. 
Apart from some local modification due to beam loading effects, as seen in (b) this wakefield 
remains robust even as the laser beam propagates though the plasma a distance of 7.5mm (as seen 
in (c) and (d)) or 5 Rayleigh lenghts. After the laser beam has propagated 2mm (as seen in (b)) 
into the plasma, one can clearly see self-trapped electrons in the first accelerating bucket. The 
radial and longitudinal localization of the self-trapped bunch is evident in part (c). After 7.5mm 
the acceleration process terminates as the depleted laser pulse starts diffracting. 
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FIG. 2: (a) A lineout of the wakefield along the z axis after 0.3mm shows that within the first bucket 
the slope of the wakefield is nearly constant and equal to eE z /(mcujp) ~ £/2 where £ = (k p (ct — z)). 
After 5.7mm of propagation (b) the wakefield has been modified by beam loading (flattening of 
the wake between 400 — 450c/u;p). This is corroborated by the p z vx z plot that is superimposed 
on the lineout of the wakefield. Pictures (a) and (b) reveal that the acceleration mechanism is 
extremely stable during the simulation. The energy spectrum after 7.5mm (c) exhibits an isolated 
spike of 0.3nC at l.bGeV with energy spread A7/7 = 3.8% corresponding to the first bucket and 
a second spike of 50pC at 700MeF with energy spread A7/7 = 1.5% corresponding to the second 
bucket. 
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FIG. 3: The normalized emittance [s]y]i = iry (Apf)(Axf) — (ApiAxi) 2 (where Api is normalized 
as indicated by the figure and the emittance is in units of Axi) and is the approximate area in 
phase space piXi. For the pictures above which correspond to the first bunch this formula yields: 
[en]x — 35TT-mm-mrad and [£N]y — 297r-mm-mrad. An upper limit for the emittance can be found 
by multiplying the typical divergences shown in the figure; this method leads to an overestimation 
which for this case is about 25%. For the second bunch of accelerated electrons (not shown in this 
figure) the emittances are significantly lower: [en]x — Wtt ■mm-mrad and [£N]y — Utt ■mm-mrad. 
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FIG. 4: _E[GeV] vs. power P[TVF] and density n/n c from equation [6]: The blue lines of constant 
power show the strong dependence of the energy of the self-trapped electrons to the density. The 



black points correspond to: (a) experiment Q], (hi experiment |n] which uses a channel lor guiding, 
(c) experiment and 3D PIC simulation {|, (d) 3D PIC simulation (| which uses a channel for 
guiding, (e) A 3D PIC simulation in [5(, (f) 3D PIC simulation presented in this article. Each of 



these points is very close to one of the blue lines indicating agreement with our scaling law. 
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FIG. 5: Picture (a) shows the laser electric field with orange on top of the electron density with 
blue after 0.4mm of propagation. The laser front is chosen sharper than the back because it has 
been found through simulations that this leads to more stable propagation. The blowout region is 
well formed. Picture (b) shows that particles at the rear of the ion channel have already reached 
velocities higher than the laser velocity and therefore are trapped. The region in z direction plotted 
in picture (b) corresponds to the z region between the broken lines in picture (a). 
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